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Abstract 

We review two subjects in the theoretical study of jet production at the Tevatron 
collider: the uncertainties in the determination of the partonic densities inside the proton 
and the uncertainties in the calculation of higher-order corrections to the QCD matrix 
elements. 
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1 Introduction 



The accurate measurements of inclusive jet distributions recently reported by the experiments 
at the Fermilab Tevatron collider [jl], |^ pose a serious challenge to the theorists' capability to 
evaluate the relative cross sections with matching precision. The current data, and in particular 
the discrepancy between data and theory at high Ej. reported by CDF , call for an improved 
assessment of the theoretical uncertainties present in the evaluation of the 1-jet inclusive 
distributions. Several groups have undertaken this task, contributing in the past several months 
to a significant progress in the field. The duty assigned to me by the organizers of the Workshop 
was to review this progress. 

We can divide the problem into two main subjects: (i) the study of the uncertainties in 
the determination of the partonic densities inside the proton (discussed here in Section 2), and 
(a) the study of the uncertainties in the perturbative calculation of the QCD matrix elements 
(discussed in Section 3). I will skip all preliminaries on jet definitions and technical aspects of 
the evaluation of the NLO matrix elements, assuming that the reader is familiar with them. 
For a complete discussion of these points the reader is referred to the references quoted below. 

2 Parton Densities 

The cross section for 1-jet inclusive production is known today up to NLO accuracy 0, ^, |^. 
It is very unlikely that full NNLO results will become available in the short term, due to 
the complexity of the calculations. A standard way to estimate the possible size of higher- 
order corrections is to evaluate the dependence of the NLO result w.r.t. variations of the 
factorization and renormalization scales. The analysis in the range 50 < ^'^(GeV) < 400 shows 
a scale dependence of the order of 10-20% |^. This uncertainty is similar to the current 
experimental uncertainties, and its small size makes the comparison between theory and data a 
compelling test of perturbative QCD, as well as an important probe of the presence of possible 
new physics. 

However accurate the perturbative evaluation of the partonic matrix elements, the ultimate 
precision of the theoretical calculations for the jet cross sections is limited by the uncertainty 
in the knowledge of the parton-density functions (PDFs) of the proton. PDF parametrizations 
are extracted from global fits to a large variety of data for which NLO QCD predictions are 
available [^. The most significant inputs are given by deep inelastic scattering (DIS), Drell- 
Yan (DY) and prompt-photon fixed-target production data. DIS provides a direct probe of the 
quark densities via the measurement of E2{x,Q'^). It also gives a value of as via the scaling 
violations of F2 at intermediate values of x (x ~ 0.3-0.5). In this x range d\ogF2/d\ogQ'^ 
receives no significant contribution from the gluon density and is still sufficiently large for 
higher-twist effects to be under control. The contributions of the different quark flavours can 
be separated by using a mixture of DIS and DY data for different beams and targets. The 
gluon density at small x can be extracted from the scaling violations of F2{x, Q"^) (for x values 
sufficiently small dF2/d\ogQ^ oc as{Q'^)G{x,Q'^)). G{x,Q'^) is rather well known today in the 
region 10~^ < x < 0.1. The most direct probe on the gluon density at larger values of x so far 
has been the measurement of large- ii^y prompt photons produced in fixed-target experiments 
with proton beams, where the Compton channel qg 57 dominates at large x over the quark- 
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Figure 1: Comparison between NLO QCD predictions for the 1-jet inclusive 
distribution and Tevatron data (figure taken from ref. ^). The PDF sets used belong 
to the CTEQ4 family J^. 



sea-suppressed annihilation channel qq — > g'j. 

The use of the NLO QCD matrix elements, together with one of the most recent global PDF 
fits (CTEQ4 P), leads to the comparison of the theoretical prediction with the Tevatron data 
shown in fig. |l[ The surplus rate found in the CDF data (not confirmed by, but nevertheless 
consistent within uncertainties with the DO measurement 0) significantly exceeds the expected 
theoretical uncertainty, estimated with the study of the scale dependence. Similar conclusions 
are obtained using independent PDF fits, such as those of the MRS group pUj] , 

Several studies have been carried out to establish to which extent the available data used 
for the extraction of the PDFs set compelling constraints on the highest-i?T tail of the jet 
distribution. Before illustrating the results of these studies, I will present a few remarks. 

The values of x that are probed and the flavour content of the initial state in jet production 
depend on the jet Ej,. In fig. I show the fractional contribution to dcr/dEj, given by the three 
possible intiial-state channels {gg, qg + qg and qq + qq). The qq channel has been separated 
into the contributions of processes which only admit s-channel exchange (namely qq — > gg and 
QiQi labelled as "QQ annih", and those which also admit a t-channel exchange (e.g. 

qiqi — ^ qifji and qiqj qiqj), labelled as "QQ t-channel". These last processes are seen from the 
figure to dominate the production at the largest Et values observed at the Tevatron, namely 
Et = 450 GeV. 

The ranges of x that are relevant to the production of jets at three different energies are 
shown in fig. ^b. While the distributions as a function of x peak at the obvious value x = Xt ^ 
2Ej,/y/S, long tails are present even at x values significantly larger than Xt- As a result, in 
order to have accurate predictions for jet production at the highest transverse energies observed 
at the Tevatron, Ej, ~ 450 GeV, the PDFs should be known for values of x up to about 0.6. 
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Figure 2: Left: relative contribution of the different initial states for the production 
of a jet at rj = as a function of the jet Ej. . Right: distribution in x for the production 
of jets with Et=200, 300 and 450 GeV. 



Since the dominant contribution to the jet rate at 450 GeV is given by quark initial states, 
it is natural to explore the possibility that the quark densities are not known with sufficient 
accuracy. First of all, given that the quark densities at large x are extracted from DIS data at 
low Q^, it is important to establish what is the range in x at low of relevance for the jet rates 
at high Q^. In fact, given the QCD perturbative evolution, it is natural to expect that most of 
the quarks found at large with x of the order of 0.5 will actually come from quarks having 
a much larger value of x at the lowest scales explored by the DIS experiments. Figure ^ 
gives the integral of the quark density above a given value of x at Q = 3 GeV (a typical scale 
for DIS data) and at Q = 225 GeV (the scale for production of a 450 GeV jet, when using the 
standard choice fi^ = ftp = Et/I). 

For X values in the range 0.5-0.6, the evolution reduces by a factor of 3 the probability 
to find a quark. This in itself is a dramatic change, which prompts some observations. To 
start with, it is important to point out that the size of the depletion at large x is controlled by 
the value of a^, since at large x d\ogF2/ dlogQ"^ oc —asiQ'^)- The larger the value of a^, the 
stronger the depletion. Therefore, although one might expect that a larger value of as would 
generate an increase in the rate of jet production, the opposite is actually true at large Et- 
This is shown, for example, in fig. |a, which is taken from the work of ref. [|r^. The fit labelled 
J (dashed line) is obtained by including the CDF jet data in the global PDF fit. In order to 
accommodate the shape of the lower-energy jet data, the fit selects a value of as slightly larger 
than that preferred by the pure DIS analysis (as(M|) = 0.120 vs. as(M|) = 0.113)0. The 
larger value of as is refiected in the lower rate predicted by this set of PDFs at high Ej.. As an 
additional comment, it should be remarked that such a significant change in the quark density, 
induced by the evolution, might be sensitive to the presence of large higher-order corrections to 
the evolution equations. In the usual analyses, the evolution is performed at NLO, consistently 

•^It is important to notice that this result is consistent with an independent evaluation of a, from the CDF 
jet data 
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Figure 3: Evolution of the integrated quark density, summed over all q and q flavours. 



with the level of accuracy of the matrix element evaluation. In the region of large x and low Q^, 
however, big corrections could come from the resummation of threshold effects, as suggested 
in [|13]. The inclusion of these effects has never been carried out, and it is hard to evaluate at 
this stage what their impact could be. 

Another interesting aspect of the evolution displayed in fig. ^ is that the number of quarks 
with X > 0.6 at Q = 225 GeV is approximately equal to the number of quarks with x > 0.7 
at Q = 3 GeV. Since the quark density drops very quickly (almost one order of magnitude in 
going from x = 0.7 to x = 0.8), we can conclude that quarks with x > 0.7 at Q = 3 GeV have 
no significant impact on the density of quarks in the region x < 0.6 at Q ~ 200 GeV. Therefore 
provided the low-energy DIS data are solid in the range x < 0.7, we expect no significant 
uncertainty due to the knowledge of the quark PDFs. The data are indeed quite accurate in 
this region, as shown in fig. ^b. For this reason, when the jet data are included in a global fit 
of PDFs the high-Erp points have no impact on the quark densities, which are almost entirely 
controlled by the small error of the DIS data points. 

Figure displays also the results of a fit obtained by removing all DIS data in the region 
X > 0.2 while including the CDF jet data (fit J'). The absence of independent constraints on 
the quark densities allows the jet data to drive them arbitrarily high, and an acceptable fit of 
the jet data is possible. However, when this fit is compared with the existing DIS data, the 
disagreement is absolutely unacceptable. The prediction for F2 obtained by the J' fit is shown 
in fig. |6, compared with the data and with the results of standard fits. At x = 0.55, F2 is 
increased by about 20-30%, which is what one would in fact expect in order to accommodate 
the 50% excess rate observed by CDF. The authors of ref. |Tl| therefore conclude that it is not 
possible to fudge the quark densities in order to fit the CDF jet data without an unacceptable 
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Figure 4: Left: Comparison between NLO QCD predictions for the 1-jet inclusive 
Et distribution and CDF data (figure taken from ref. J71[/J. The dashed line uses a 
PDF fit obtained including the CDF jet data in the fit. The dash-dotted line uses a 
PDF fit which includes the CDF data but leaves out all DIS data for x > 0.2. Right: 
Comparison between the PDF fits used in the left figure and BCDMS data for F2 (figure 
taken from ref. ^Jl|/j. 
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Figure 5: Left: ratio of the CTEQHJ and CTEQ4M gluon densities. Right: com- 
parison between NLO QCD predictions for the 1-jet inclusive Ej. distribution using the 
CTEQHJ set and the Tevatron data (figure taken from ref. ^). 



clash with the high-precision DIS data. 

An alternative approach was pursued by the CTEQ group. These authors exploited the 
fact that the experimental and theoretical uncertainties present in the analysis of the prompt 
photon data might leave enough room to accommodate larger gluon densities at large x. From 
fig. 1^ we see that in order for processes with initial-state gluons to explain the 50% excess 
rate, the gluon density should increase by at least a factor of 2. The CTEQ group performed 
a global fit [p!4| , ^ including the high-i^^ CDF data and assigning to them a larger weight in 
the determination. The gluon density obtained in this way, called CTEQHJ, is compared 
with a "standard" gluon density in fig. ^a. The global of the fit is marginally worse than 
the best (a detailed discussion of this point can be found in ref. In particular, the 

theoretical predictions for the prompt-photon spectrum obtained using this set of PDFs are 



perfectly consistent, within uncertainties, with the experimental data |T^. The comparison 
between QCD and data for the iet-E^ distributions is shown in fig. ^b. As can be seen, most 
of the excess rate has been removed by this choice of PDFs. 

Before concluding this section I would like to make a few more comments. During the 



summer, new results on neutrino DIS from the CCFR collaboration were released 151. The most 



significant result of these analyses is a new extraction of as from the large-x scaling violations 



in F2 and xF^. Contrary to all previous extractions of from DIS 1 16], this measurement 



is in perfect agreement with the values of as extracted at higher energies (e.g. at LEP) or 
in r decays |T^. In particular, the new CCFR measurement of as{Mz) = 0.119 ± 0.005 is 



significantly larger than the previous average of DIS measurements, as{Mz) = 0.112±0.004 [18 
It is hard to evaluate what the impact of these data will be on the calculation of jet cross sections 
for the Tevatron. The F2 and xF^ data have not been released as yet, so that no global fit of 
PDFs and their evolution is available. It is extremely encouraging that the new DIS as value 
is consistent with that extracted from the analysis of the jet data |T^. Even more so because 
these determinations are now consistent with the extraction of as from other independent 



measurements |T^. It is likely that, because of the enhanced evolution of the large-x valence 
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densities, the high-i^y CDF anomaly will be even more significant, as we already discussed 
above. As the CTEQ studies have shown, however, this could be absorbed into an acceptable 
change in the gluon densities. Independent measurements of the large-x gluon component of the 
proton are at this point of extreme importance. Higher statistics prompt-photon production 
data are becoming available [BO], and will hopefully help to settle this issue soon. 



3 Higher orders in perturbation theory 

As mentioned in the previous section, the study of the scale dependence in the range 50 < 
Ej.{GeV) < 400 shows an effect of the order of 10-20% 0, It is reasonable to expect that 
the NNLO corrections, in this Ej, range, will be of this size. For larger Ej, values the scale 
dependence becomes more and more pronounced 0, ^j, indicating that NNLO effects can no 
longer be neglected. Whether the actual threshold for the onset of significant NNLO contribu- 
tions is 400, 450 or 600 GeV, is clearly a question of extreme relevance for the interpretation 
of the CDF data. 

To answer this question we must find a more solid way of estimating the higher-order 
corrections than just looking at the scale dependence of the lower-order results. Direct estimates 
of higher-order corrections are possible when these are dominated by terms whose structure does 
not depend on the specific details of the hard process under consideration. A typical example 
is given by the large logarithms that are associated to the emission of collinear gluons from 
the initial state of a hard process. These logarithms are independent of the details of the hard 
process, and can be resummed to all orders of perturbation theory, giving rise to the standard 
DGLAP equation which governs the evolution of the parton densities inside the proton. 

The resummation of higher-order contributions is generally possible when two different hard 
scales (say fi and Q, with Q ji) are present in the problem. Gluons whose wavelength is 
of the order oi 1/ ^ cannot be sensitive to phenomena taking place at distances of the order 
of 1/Q. If their emission gives rise to large corrections of order as\o^{Q/ ^) {k < 2), as is 
often the case, these corrections are therefore universal and can be resummed (for a thorough 
discussion of these issues at NLO, and for a complete set of references to the relevant literature, 
see refs. |^). In the case of high-energy jet production, the two different hard scales 
correspond to the energy of the jet itself and the maximum amount of energy that can be 
released by initial-state radiation. The closer one gets to the threshold Mjj Shad (where 
Mjj is the invariant mass of the dijet system and y/ Shad is the total energy in the pp CoM 
frame), the smaller the amount of energy available for radiation. Given that parton densities 
fall steeply, the phase space available for soft radiation is already significantly reduced before 
T = Mjj/ Shad approaches 1. In the case of DY production in fixed-target experiments, for 
example, large corrections are already found for values of r of the order of 0.5 [E^, It is 



believed that higher powers of aslog^(l — r) are responsible for the increased sensitivity of 
the jet cross section at high Ej., and that these are the leading corrections that need to be 
resummed to provide a better estimate of the jet cross section beyond NLO. 

A technical discussion of how the resummation of threshold corrections is performed, in 
addition to references to the original literature, can be found in ref. |T^. Here I will limit 
myself to presenting one result of interest for the study of high-energy jet production. Rather 
than considering the resummation corrections to the Ej, distribution of jets, I will consider the 
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case of high-mass dijet pairs. The reason is that resummation is simpler in this case, since for 
this variable all Sudakov effects are present only in the evolution of the initial state, so that 
the formalism is almost identical to that of the DY processes. The study of the resummation 
for the 1-jet inclusive Ej. distribution is in progress. 
In fig. M I show the following quantities 



5(3) 5(3) x(3_) 5(3) , 5(3) , ^(3) 

"gg "qg "qq "gg ^ "qg ^ "qq 
a(2) ' a(2) ' a(2) ' a(2) 



(1) 



where S^^^ is equal to the resummed hadronic cross section in which the terms of order have 
been subtracted, and cr*^^^ = + a^J + a^^ is the full hadronic LO cross section (of order 
ag^). I use as a reference renormalization and factorization scale /i = Mjj/2. Notice that for 
large invariant masses the effects of higher orders are large. To understand how much is due to 
the first-order corrections (which are exactly calculable p3| , and how much to corrections of 
order and higher, I show the following quantities in fig. |^: 

5(4) 5(4) .(4_) 5(4) , ^(4) , ^(4_) 

"gg "qg "qq "gg ^ "qg ^ "qq /^^n 
(t(3) ' (t(3) ' a(3) ' (t(3) ' 

where 5^^^ is now equal to the resummed hadronic cross section with terms of order sub- 
tracted, and cr''^^ is an approximation to the full NL cross section, summed over all subprocesses, 
obtained by truncating the resummation formula at order ag^. This figure shows that most of 
the large K factor is due to the pure NLO corrections, with the resummation of higher-order 
soft-gluon effects contributing only an additional 10% at dijet masses of the order of 1 TeV. Only 
above Mjj > 1200 GeV, corresponding approximately to > 600 GeV, are the resummation 
effects non-negligible. 

These results should only be taken as an indication of the order of magnitude of the correc- 
tion, since we have not included here a study of the resummation effects on the determination 
of the parton densities. As was already mentioned previously, it is quite possible that resum- 
mation effects significantly influence the determination of the large-x structure functions from 
low energy data. Equally important effects could appear in the evolution from low to high Q^. 
It would therefore be important to reexamine the extraction of the large-x non-singlet structure 
functions, in the light of the resummation results for the DIS process, before firmer conclusions 
can be drawn on the significance of the present jet cross section discrepancy. 

From this preliminary study it seems, however, unlikely that the full 30-50% excess reported 
by CDF for jet E^s in the range 300-450 GeV could be explained by Sudakov resummation 
effects in the hard process. 

Before concluding this section, I want to mention a puzzling observation recently made in 



ref. ||26|. The authors studied the scheme dependence of the jet cross section, by using the MS 
and DIS versions of the same PDF fit, folded with the NLO jet cross section evaluated in the 
two respective schemes. The two evaluations should yield equal results, up to terms of NNLO. 
However, the differences turn out to be numerically large, of the order of 40% at the highest 
Et's. The DIS calculation, in particular, is in good agreement with the CDF data. Whether 
this difference corresponds to a genuine scheme dependence, and therefore a true uncertainty 
due to the ignorance on corrections beyond NLO, or whether it is an artefact of the way the 
DIS PDFs were extracted or evolved, this is still not clear. 
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Figure 6: Contribution of gluon resummation at order a^^ and higher, relative to the 
LO result, for the invariant mass distribution of jet pairs at the Tevatron. 
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Figure 7: Contribution of gluon resummation at order a s'^ and higher, relative to the 
truncated O {a s^) result, for the invariant mass distribution of jet pairs at the Tevatron. 
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4 Conclusions 



The main conclusions of this review can be summarized as follows: 

• The NLO approximation is sohd up to Et — 500 GeV, with uncertainties at the level 
of 10-20%. The effects of higher-order corrections, estimated here in the case of dijet 
production by resumming the largest leading soft logarithms appearing as x — > 1, will 
become non-negligible only for transverse energies larger than 600 GeV. 

• Modifications, at large x, of the gluon density evaluated so far seem possible without 
affecting the overall quality of global PDF fits. It seems possible to absorb a large fraction 
of the excess of high-i?y jets reported by CDF, by allowing for a harder gluon at large x, 
without affecting the comparison of theory with other sets of data. 

• While resummation corrections do not seem important in the evaluation of the hard cross 
sections, their effect could be non-negligible on the evolution of the valence quark densities 
from the low-Q^ region (where they are measured) to the high-Q^ region relevant for the 
Tevatron jets. 

• The new measurements of F2 and xF^ by CCFR call for new analyses of the global fits 
to the PDFs. Only at that point will a new assessment of the situation with the high-E'r 
Tevatron jets be possible. 

The last two points indicate that, in spite of the significant progress we witnessed recently, 
there is still room for new ideas and for improvement. 
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